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Abstract. The Ua{1) problem of QCD is inevitably tied to the infrared behaviour of quarks and gluons with its 
most visible effect being the 77' mass. A dimensional argument of Kogut and Susskind showed that the mixing of the 
pseudoscalar flavour-singlet mesons with gluons can provide a screening of the Goldstone pole in this channel if the 
full quark-quark interaction is strongly infrared singular as ~ 1/fc*. We investigate this idea using previously obtained 
results for the Landau gauge ghost and gluon propagator, together with recent determinations for the singular behaviour 
of the quark-gluon vertex. We find that, even with an infrared vanishing gluon propagator, the singular structure of the 
quark-gluon vertex for certain kinematics is apposite for yielding a non-zero screening mass. 
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1 Generalities on the ^7^(1) Anomaly 

It has now been long accepted that the strong interaction is 
described by a Yang-Mills theory, whose non-Abelian charac- 
ter exhibits the properties of asymptotic freedom and confine- 
l ment. However, despite our accedence to this we still have no 
satisfactory understanding of the confinement phenomena it- 
self. It is envisaged that it is in the infrared structure of QCD 
where this mechanism lies, and so we must probe in a momen- 
tum region characterised by large values of the coupling. Here, 
standard tools such as perturbation theory are inapplicable and 
lattice calculations typically use volumes too small to reliably 
probe the deep infrared. 

The idea of 'infrared slavery' gave rise to two possible be- 
haviours for the strong running coupling at small momenta 1 1 1 : 
either it exhibits an IR fixed point or diverges at the origin. 
From functional methods in Landau gauge it is the former be- 
haviour that seems likely in the Yang-Mills sector |!2l[3]|4l|5]|6l . 
However, such a fixed point behaviour of the coupling does not 
predicate that the vertex functions themselves are infrared fi- 
nite. Indeed, depending upon the number of external ghost and 
gluon legs, many IPI Green's functions are found to be sin- 
gular in the infrared |l5]|6]|2|. For the quark-gluon vertex, this 
issue is addressed in H). In this paper it is our aim to address 
yet another puzzle of QCD, with the infrared behaviour of the 
quark-gluon vertex playing a crucial role. 

We know that the up and down quark are very light, and so 
may consider them massless to good approximation. To a lesser 
extent, the same may be assumed for the strange quark and so 
we expect the Lagrangian of QCD to exhibit an approximate 
SU{Nf) X SU{Nf) symmetry, with Nf = 3. Assuming the 



chiral limit is realised, this is broken dynamically to a diagonal 
SUv{i), and so gives rise to eight Goldstone bosons. 

However, we also expect the same mechanism to break the 
C/y!i(l) symmetry of the QCD Lagrangian. Thus one should 
find a ninth isosinglet pseudoscalar boson corresponding to the 
spontaneous breaking of this symmetry. However, on exam- 
ining the spectrum of observed mesons this proves to be one 
Nambu-Goldstone boson too many. The lightest candidate for 
this is the rf, whose mass of m^/ ^ 958 MeV is far from being 
'light'. In fact, it has been shown that if the C//i(l) symmetry 
is not explicitly broken, then the mass of the rj' must be less 
than a/S^tt ^ 250 MeV f9]. Even allowing for the mass of 
the strange quark, the 77' remains far too heavy to be solely de- 
scribed in this manner. This is the U{1) problem of QCD. 

A step towards resolution came from the recognition that 
the classical C/a(1) symmetry is anomalous, i.e. broken by quan- 
tum mechanical effects. However, it is found that the associated 
current is a total divergence, and hence no symmetry break- 
ing contribution is obtained for finite order in perturbation the- 
ory |10|. Thus, the mechanism of the anomalous breaking of 
the ?7yi(l) symmetry must be wholly non-perturbative in na- 
ture. 

There have been several suggestions as to how the rj' ob- 
tains its unexpectedly large mass. In the early days of QCD 
Kogut and Susskind pointed out, on dimensional grounds, that 
a contribution to the rj' in the chiral limit could be obtained by 
the mixing between two infrared enhanced gluons, with mo- 
mentum space propagator D{k) ^ \/k'^ for ^ Q llll . A 
few years later topological solutions termed instantons were 
discovered in QCD [12|. In contradistinction to the confine- 
ment driven Kogut-Susskind scenario, these gave rise to an al- 



2 



Alkofer et al: Ua(^) anomaly and rj' mass from an infrared singular quark-gluon vertex 



ternative solution of the Ua{^) problem: Instantons were shown 
by 't Hooft to lead to the non-conservation of the axial charge 
and so induce a 2Nf fermion operator giving rise to a non- 
zero Ty'-mass in the chiral Umit Ml 311 14 1 . Again a few years later 
Witten and Veneziano II15II16| [T71| proposed their solution of 
the problem by considering an expansion of QCD in Nf/Nc, 
where Nf and Nc are the number of flavours and colours re- 
spectively. They showed that the correct pattern of the C/^(l)- 
symmetry breaking can be obtained if the anomalous mass of 
the rj' is related to the topological susceptibility of the the- 
ory. Since this susceptibility is not necessarily generated by 
instantons Veneziano named his paper 'C/(l) without instan- 
tons' |[T6 |. Indeed, a different type of topologically non-trivial 
gluon field configurations leading to a non-vanishing topolog- 
ical susceptibility are center vortices as e.g. discussed in ref. 
108 1 . 

Certainly, with the Witten- Veneziano mechanism at hand 
confinement and topology based explanations for the Ua{^) 
anomaly are not mutually exclusive. The reason is simply that 
well-established suggestions for topology driven mechanisms 
for confinement exist |50|. Thus the same topologically non- 
trivial gauge field configurations could be responsible for both 
confinement and the resolution of the L/yi(l)-problem. 

In this letter we concern ourselves with the Kogut-Susskind 
mechanism employing an approach similar to that given in |fT9l 
l20l . We employ solutions to the Dyson-Schwinger equations 
for the propagators of QCD and the quark-gluon vertex to de- 
termine the anomalous mass of the rj' in the chiral limit. The 
technical details are discussed in sections |2] - 14] We present 
our result in section |5] and discuss its relation to the Witten- 
Veneziano mechanism in section |6] We summarise in section 
|7] We wish to note that all calculations are done in Euclidean 
momentum space. Furthermore, we will ignore all effects from 
isospin breaking and work in the isospin limit. 



2 Flavour mixing 

A suitable basis to describe the flavour content of mesons is 
given by the SU (3) singlet and octet basis: 

|7r°) = {\uu) - \dd)) /V2 

|778) - {\uu) + \dd) - 2 /%/6 (1) 

|?/o) = {\uu) + \dd) + \ss)) /\/3 

The octet-singlet mass squared matrix in the isospin limit is 
thus given by 







with matrix elements 

M2 = 
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Here, we have employed m^s = 2m\ — m^ from the Gell- 
Mann-Oakes-Renner relation to make the substitution on the 
right-hand side. With the basis of ([T]| the pion is decoupled in 
the isospin limit. Thus we concentrate on the 2x2 sub-matrix 
that mixes the rj and rj' 1161 : 



imi 



2V2{ml 
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(4) 



The anomalous mass term, is related to the topological 
susceptibility through the Witten- Veneziano formula IITSI 

AT . 

(5) 



9 9 



which includes the pion decay constant /q ~ 93MeV. By now 
its derivation has been refined and it has been adapted to lat- 
tice gauge theory and found to nicely agree with theoretical ex- 
pectations and experimental data |21 ,22,23,24,25], with val- 
ues reported of the order ^ (191 ± SMeV)''. It is impor- 
tant to note, however, that the corresponding Witten- Veneziano 
mechanism does not make any reference as to which degrees 
of freedom are responsible for the generation of x^- We come 
back to this point in section |6] 

The 770 and 778 are not the physical mass eigenstates, and 
so we diagonalise the 2x2 mass-matrix by introducing the 
following unitary transformation: 



(77 ?/) 



cos ( 
sin^ 



— smf 
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(6) 



with d the singlet-octet mixing angle. This gives rise to the fol- 
lowing mass eigenstates: 
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(7) 



with m„ corresponding to the 77' and 77 respectively, and 
the shift A given by: 



A^d{ml 



mlY -"^{ml- ml) + ^ 



771' 
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A phenomenological value for the mixing angle can be calcu- 
lated from: 



tan 9 
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(9) 



with the current values mx ~ 498 MeV, m.,^ = 135 MeV | 
Here, the required input is the anomalous mass contribution 
777^. The Random Instanton Liquid Model provides a mixing 
angle of 6* ~ -11.5° II27II28L whilst knowledge of the 77 and 77' 
masses, respectively about 547 and 958MeV, suggests ^ —20°. 
A detailed study of the available observables yields an average 
value of -15.4° |29|. 



3 Anomalous Contribution of the Gluon 

If we ignore the contribution from the gluon anomaly, that is by 
setting m\ ~ 0, we would expect the flavour content of our 77, 
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rf to be analogous to that of the uj, p system where ideal mix- 
ing is realised to very good approximation. Diagonalising ^ 
entails that the mass of the eta be degenerate with that of the 
pion, with flavour content (mm + dd) / y/2, whilst the 77' would 
be a pure ss pseudoscalar. We know this to be far from the 
truth since this would provide a mixing of 6* ~ —54.7°, hence 
the [/4(l)-problem and the desire to include anomalous contri- 
butions from the gluon. 

To do so in the Dyson-Schwinger and Bethe-Salpeter ap- 
proach is a formidable task, for presently we are limited by ker- 
nels that include only re-summations of gluon ladders. These 
do not include the necessary annihilation (and consequent re- 
combination) of qq pseudoscalar meson into two gluons. Typ- 
ically, the C/^(l)-breaking term is put in by hand and adjusted 
as a free parameter to obtain phenomenological results. 

The minimal diagram that includes the gluon Adler-Bell- 
Jackiw (ABJ) diagram is given in fig. [T] We later comment 
upon the contribution of diagrams with more than two gluon 
exchange. 



k - pp. 




k + P/2 



Fig. 1. The diamond diagram n{P^). The crossed gluon exchange 
contribution gives rise to an additional factor of two. 



For a simple model, Kogut and Susskind gave arguments as 
to how a non-vanishing contribution to the ry'-mass in the chi- 
ral limit could be obtained 1 1 1 1 . This has been explored phe- 
nomenologically in a previous study 1 19,20|, in order to gauge 
the magnitude of the effect should it be realised. There, an 
ansatz for the gluon propagator has been employed which is 
infrared singular supplemented by a perturbative tail. Though 
such a form cannot naturally arise from consideration of gluons 
alone, at the time there were already suggestions that in Lan- 
dau gauge the ghost sector was strongly infrared enhanced and 
could thus enter into the effective gluon propagator implicitly. 
Using parameters given in that paper, a screening-mass squared 
of about half the required amount is obtained, giving credence 
to the mechanism. However, nothing definitive may be claimed 
until dressed Green's functions are employed, whose behaviour 
is well-founded in the infrared. 

In order to move beyond this qualitative study, we use ex- 
plicit calculations for the gluon obtained in previous calcula- 
tions |31 32|. Since these prove to be infrared vanishing, it is 
important that we have a self-consistent and motivated dressing 
for the quark-gluon vertex that compensates for this precisely, 
again leading to a non-vanishing topological susceptibility in 
the chiral limit. Our starting point for the calculation is the di- 
agram shown in fig. [T] It consists of two anomalous AW tri- 



angle diagrams 113311 . which here are obviously non-Abelian in 
nature fT3l. These triangles are linked via the exchange of two 
gluons. Note that should we exchange three gluons or more in 
the i— channel, one would couple in a pseudoscalar J^'~"' — 
^ glueball, whose large mass of 2 GeV Ii34l would heav- 
ily suppress the contribution. For multiple gluon exchange in 
the s— channel, however we essentially re-sum a gluon ladder 
which is related to meson exchange [35J. For this, the scalar 
part of the vertex which we have neglected to model here is 
important. 

Staying with two-gluon exchange and thus restricting our- 
selves to a qualitative study, we write the amplitude of fig. [T] 
explicitly as an integral over the gluon momenta: 

7T(P2) = I TkGl^^{P,k)DZ{k+)DZ{k-)Gt{~P,k) , 

with dk shorthand for the invariant measure d'^k/{2Tr)*, and 
k± = k ± P/2. We have chosen an equal momentum partition 
to simplify the resulting equations; the final result should be 
independent of this choice. The quark triangle, G^fl, arising 
from the axial anomaly is: 



G';!i{p,k) 



dqi^T 



r,,{P,q)S{q+)n{q+,q-k) 



xSiq-k)r^{q^,q-k)S{q-) 
+tr r„iP,q)S{q+)r^{q+,q + k) 

xS{q + k)r^{q-,q + k)S{q.) 



(11) 

where the trace is over Dirac and colour indices, and the sec- 
ond term accounts for the crossed gluon-exchange. Solution of 
this equation requires knowledge of the pseudoscalar Bethe- 
Salpeter amplitude (BSA) for the Pr^iP, k), the full quark prop- 
agator S{k), and the quark-gluon vertex r^{k, q). 

4 Calculational Inputs 

To calculate the quantities required for computation of the di- 
amond diagram, we employ the Dyson-Schwinger approach to 
QCD, see 1 36, 37,38,39] for reviews. These comprise of an in- 
finite tower of coupled integral equations that interrelate the 
fundamental Green's functions of the theory. Despite analytic 
studies being possible in certain kinematical situations |[3]|5]|6] 
40, 7 1, in principle one must employ a truncation scheme and 
perform the calculations numerically. In the present context it 
is, however, important to mention that anomalous diagrams are 
represented exactly if the underlying symmetries are respected, 
see e.g. refs. L41..42J . 



4.1 Yang-Mills Propagators 

The pure Yang-Mills (YM) sector of QCD is comprised of both 
ghosts and gluons, whose propagators in Euclidean space are 
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written: 



Similarly, the gluon dressing function may be fitted ||32l : 



D 



ah 



sab 



Pf^pA Z{p 



(12) 

where G{p^) and Z{p'^) are the ghost and gluon dressing func- 
tions respectively. These satisfy the Dyson-Schwinger equa- 
tions shown pictorially in fig. |2] and have been studied in Lan- 
dau gauge by employing truncations at the level of the vertices. 
It is found that in the infrared these can be described by power 
laws Il2ll3,44|: 



G{p')-{p'Y 



Z{p') - {p'f 



(13) 



with K a positive constant, giving rise to an infrared diverging 
ghost propagator and vanishing gluon. Moreover, with a bare 
ghost-gluon vertex one finds k = (93 - V1201) /98 ~ 0.595, 
and so we see IR dominance of the ghost in Landau gauge. This 
fact greatly affects our calculation of ( fTOt since we now entirely 
depend upon the dressing of the quark-gluon vertex to yield the 
sufficient infrared enhancement that will lead to non-vanishing 
7T(P2 _> 0). 

From the ghost-gluon vertex we may define the running 
coupling: 

aip') = a^G\p^)Z{p^) , (14) 

which may be parameterised such that the numerical results for 
Euclidean scales are accurately reproduced |32|: 



a&t{p^) 



«.(0) 



An 



P 



p'^/A^ 



QCD 
1 



In(p2/yl2 



QCD^ 



QCD +P 



(15) 



Here /3o = (HiVc - 2A^/)/3, and as(0) is the fixed point in 
the infrared, calculated to be 8.915/A*'c for our choice of k. 




t2 J- /12 
+ ^^QCD 



(16) 



such that no further singularities are introduced. Here, 7 = 
(-13iVc + 4Nf)/{22Nc - 47V/) is the one-loop value for the 
anomalous dimension of the gluon propagator, and = 0.2 
at the renormalisation scale fi^ — 170 GeV^. We use ^qcd = 
0.5 GeV'^ similar to the scale obtained in ref. [i32i. 



4.2 The Gap Equation 

A far simpler system to that of the Yang-Mills sector, though 
one of equal if not greater importance, is that for the quark 
propagator, shown in fig. [3] This may be written explicitly: 



s-\p)^Z2[s("\p) 

Z\ Zi 



(17) 



y d^k-f^,s{k)r,{k,p)D^,{p - k) 



As is evident, this is a non-linear equation dependent upon the 
form of the quark propagator itself, the gluon propagator and a 
higher order Green's function, the quark-gluon vertex r,^{k,p). 
The quark propagator consists of a Dirac-odd and Dirac-even 
part, and so is completely described by two momentum depen- 
dent functions: 



(18) 



with A/(p^) the mass function and Zf{p'^) the wave-function 
renormalisation, which are both assumed to be non-singular in 
the infrared. This behaviour is indeed what is seen in numerical 
calculations lISTU Sl. In the following we consider the quarks in 
the chiral limit, i.e. [^^"•'(p)]^^ — i-p. Then M{p^) only con- 
tains effects from the dynamical breaking of chiral symmetry 
and in the — Ty'-mixing we only have the contribution m^j 
from the topological susceptibility 




Fig. 2. The gluon and ghost Dyson-Schwinger equation. Filled blobs 
indicate dressed propagators and vertices. 



4.3 Quark-Gluon Vertex 

Earlier we stated that in light of the vanishing nature of the 
gluon propagator in the deep IR, it is the behaviour of the 
quark-gluon vertex there that is crucial. 




Fig. 3. The quark Dyson-Schwinger equation. Filled blobs indicate 
dressed propagators and vertices. 
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4.3.1 Vertex for all scales vanishing 

From the IR analysis of Il40l l8l we know that the quark-gluon 
vertex, r'" = V-^M' should have the behaviour: 

r^ipY^-'^" , (19) 

when the sum of all incoming momenta, p, are vanishing in the 
infrared; that is, it is strongly infrared enhanced. They find that 
all twelve possible structures of the quark-gluon vertex may 
contribute in the infrared, with vector and scalar pieces playing 
a dominant role. Thus one may parameterise r"^ by: 

r„ [p, q) = Li {p, q)-iu-i L3 (p, q) [p + q)^ (20) 

with p and q the incoming and outgoing quark momenta re- 
spectively, and k — p — q the outgoing gluon momentum. 

We wish to find a non-vanishing n{P^ 0), and thus 
look at the quark-triangle G'^^'l in the limit of the incident gluon 
momenta vanishing, i.e. k,P ^ j^qcd- We note that we may 
write: 

Gf. - Is-'e^.^^P'^k" I (P2, fc^ fc . P) , (21) 

where / represents the internal loop integral of the quarks, ob- 
tained by separating out the tensor structure. We then evaluate 
the behaviour of / numerically, using the power laws (fTST i and 
( fT9] l for the gluon propagator and the quark-gluon vertex, re- 
spectively, and an infrared finite quark propagator and Bethe- 
Salpeter amplitude. We then find that the dominant infrared be- 
haviour of the triangle is given by: 

/(F2,0,0)^ . (22) 

On dimensional grounds this behaviour is insufficient to gen- 
erate 77 7^ in the limit P^ — > 0. We thus conclude that the 
uniform scaling behaviour fT% of the quark-gluon vertex is not 
sufficient to account for the Kogut-Susskind mechanism. 

4.3.2 Soft Singularity 

However, there is not just one momentum configuration giving 
rise to a singular structure in the vertex. There also exists a 
soft collinear-like divergence dependent only upon the external 
gluon momentum k^: 

P~(fc2)-"-i/' (23) 

This additional divergence has been identified from the Dyson- 
Sch winger equation for the quark-gluon vertex in ref. fE]. Its 
appearance in the quark-antiquark scattering kernel for heavy 
quarks leads to a linearly rising potential and consequently to 
quark-confinement. Since the derivation of eq. ( |23] | is some- 
what involved we do not wish to repeat it here but instead refer 
the interested reader directly to ref. |8 1. 

For the integrand 7 at small momenta this soft divergence 
leads to 

/(P2,0,0)- (P^)-""-! , (24) 



a behaviour which we have obtained both by analytical power 
counting and also numerically. To gauge the impact of this re- 
sult on our evaluation of 7T(P^ 0), we insert {2l[ into (fTOb 
to obtain: 

n{p')^ujTk{k-pf ^g±l ^1^72 . (25) 

A simple check of the dimensions indicates that a non-zero 
contribution to the anomalous mass is obtained. Only the mag- 
nitude of this needs now be determined. To this end, we must 
still determine the Bethe-Salpeter amplitude for the 77. 

4.4 The Bethe-Salpeter equation 

An essential input to the diamond diagram is the Bethe-Salpeter 
amplitude for the -q. This can be obtained by solving the pseu- 
doscalar Bethe-Salpeter equation (BSE), which describes qq 
bound-states: 

Ptuip-.P)^ lTkKl%ip,k;P)[Sik+)rik;P)Sik^)l^ , 

(26) 

where K represents the (unknown) quark-antiquark scattering 
kernel, fc+ = fc + rjP, fc_ — k + {rj — 1)P, rj is the momentum 
partitioning of the quarks. The meson has total momentum P, 
taken in the rest frame, with P^ = — and m the mass of the 
bound-state. 

Fundamental to any study of pseudoscalar mesons is the 
axial-vector Ward-Takahashi identity (AV-WTI): 

P^P5;(fc;P) = 5-i(fc+)iA^z75 + lx'}ij,S-'{k^) 

^ ^ (21) 

-Mcir^{k;P)-ir^{k;P)Mc . 

Since the left-hand side involves the BSE, whilst the right- 
hand side involves the quark SDE, a non-trivial relationship be- 
tween the kernels is established. Thus one must treat the Bethe- 
Salpeter kernel K and the quark-gluon vertex in such a way that 
the AV-WTI is satisfied, otherwise a massless pion will not be 
found in the chiral limit. A generalisation of this identity to the 
case of the flavour singlet channel has been given in ref. ||43l . A 
suitable truncation scheme used frequently |38| is the rainbow- 
ladder approximation which we also employ here. 

4.5 Phenomenologically qualitative model 

In rainbow-ladder approximation the quark-gluon vertex is re- 
stricted to its part. We therefore simplify the parametrisation 
(|20] | further and write: 

P,(fc,p) = yl2)Li(fc^^2)7, . (28) 

Here the vertex renormalisation constant Zip provides multi- 
plicative renormalisability of the vertex and also of the result- 
ing quark SDE and pseudoscalar BSE. The vertex dressing Li 
is parameterised by: 

L (z) ^ [ ^ / di zd^ 

\z + d2) \l + z/d2 dl + iz-d^f 
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z 



in 



Poctfj. \log{z/d2) z-d2 



-2S ' 



(29) 



with z — k'^ the gluon momentum and di sets the strength 
of the interaction in the infrared. The scale ^2 is set to ^2 = 



{O.bGeVy 



/12 

^^QCD' 



i.e. similar to the scale found in the 



Yang-Mills-sector. The parameter d^ determines the size of a 
term added to give additional integrated strength in the inter- 
mediate momentum regions. 

This construction is necessary to produce meaningful re- 
sults: Since the AV-WTI forces us at present to only consider 
the 7^-part of the quark-gluon interaction we necessarily miss 
some interaction strength in the intermediate momentum re- 
gion. This cannot be compensated for by an increase of di as 
discussed in the following. 



5 Results 

Naively, the diamond diagram consists of a twelve dimensional 
integral which can be reduced to eight dimensions by symme- 
try considerations. In the limit the integrand simplifies 
further and the problem may be reduced to a five dimensional 
integral. We solved this integral both, by Monte-Carlo tech- 
niques and using adaptive quadrature routines and found the 
latter method to be superior. Sufficient accuracy of the order of 
one percent is obtained within a few minutes of CPU-time. 

In table [T] we give our results for the calculation of the di- 
amond diagram shown in fig. [U employing the model of the 
quark-gluon interaction discussed above. Within the limitations 
of the model we may not expect quantitatively coiTect results. 
Nevertheless we obtain reasonable values for the masses and 
the mixing angle at least for the last parameter set considered. 
For this set, we find a topological susceptibility. 



X 



= (169 MeV)"* 



(30) 




P3' [Gev'] 



Fig. 4. Our model dressing for the 7^ part of the vertex. See the text 
for a description of each. 



In fig.|4]we plot two examples of Li. The parameters of the 
dashed line have been fitted such that the physical masses of 
the TT and the p are reproduced, employing di and ^2 only. Here 
we need a large value of di — 13, which leads to an extremely 
strong infrared singularity; employing such a form in the dia- 
mond diagram would yield results for the topological mass m\ 
which are orders of magnitudes too large to account for the ex- 
perimental 1] and yy'-masses. To account for this, we introduce 
extra interaction strength in the intermediate momentum region 
parameterised by ^3. This additional bump generates correct tt 
and p-masses whilst allowing for a reduction of the infrared 
scale di to a realistic order. Note, that this particular shape and 
strength of interaction is also suggested by the results for the 
running coupling of the quark-gluon vertex, found in [8J . 

Once our parameters are fixed, we obtain solutions for the 
Bethe-Salpeter amplitude in the chiral limit, which are then em- 
ployed in the calculation of the diamond diagram. 



in qualitative agreement with lattice results [25 1. Qualitatively 
we do see that it is precisely the infrared collinear behaviour 
of the quark-gluon vertex that is necessary to obtain a non-zero 
contribution to the topological susceptibility. 

A quantitative study would require a self-consistent solu- 
tion for the quark-gluon vertex employing fewer approxima- 
tions. Then, at least in the ^ limit, one need not solve the 
pseudoscalar BSE since the leading term ^5B{q^)/ f.^ in chiral 
limit is known from the quark DSE. We would then need to 
consider higher-order diagrams corresponding to the diamond 
diagram, with n-gluon exchange taking place in the s-channel. 
This ladder re-summation of gluons is often represented by 
a single meson exchange. This picture has been explored in 
ref. 1.351 as a means of generating the topological mass of the 
f]'. However, since the actual mechanism of anomalous C/4(l) 
symmetry breaking may not be explained in terms of gauge in- 
dependent particle exchange, such as a pion, this picture cannot 
be entirely adequate. 

6 Discussion of topological effects and 
infrared singularities 

As we have seen the infrared divergence of the quark-gluon 
vertex plays an important role in a confinement-based explana- 
tion of the rj' mass and thus the ?7a(1) anomaly. On the other 
hand, it is evident that the 77' mass is linked to the topological 
susceptibility via the Witten-Veneziano formula (|5]l and thus 
topologically non-trivial gluon field configurations. How can 
such a dichotomy be resolved? 

In a first step to a potential solution of this puzzle we note 
that it is not only instantons that are capable of providing a 
non-vanishing topological susceptibility. As a matter of fact, 
there is evidence [45,46 47 1 that on the lattice, field configura- 
tions with non-trivial holonomy have been misidentified as in- 
stantons. They are more like Kraan-van-Baal-Lee-Lucalorons 
[|48,49[ than instantons. And, amongst these and other topolog- 
ically non-trivial gluon field configurations, also center vortices 
(see e.g. the reviews [50]) provide a non-vanishing topological 
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Table 1. Numerical results of the calculation of the diamond diagram for a range of parameters. The n and p phenomenology is only marginally 
affected by these small changes in infrared strength, yet note the strong effect on the calculated anomalous mass. 



susceptibility II18II511 . At least, for these latter field configu- 
rations we know that they provide an explanation for quark 
confinement II50II52I . and that in most gauges, including the 
Landau gauge, a substantial amount of these configurations 
live on the Gribov horizon. The latter property is responsible 
that, when removing center vortices from a lattice ensemble, 
not only the string tension vanishes but also the Landau gauge 
ghost propagator becomes infrared finite [53 1. 

In the confinement based scenario described within this let- 
ter, the infrared divergence of the quark-gluon vertex is the 
cause for the 77' mass. It is driven by the infrared divergence 
of the three-gluon gluon vertex and thus eventually by the in- 
frared behaviour of the ghost propagator |40 |. In Landau gauge 
(as well as in some other gauges) this infrared divergence is, 
according to the Gribov-Zwanziger scenario, assumed to be 
caused by the dominance of field configurations on the Gri- 
bov horizon (or, more precisely, on the non-vanishing over- 
lap of the boundary of the fundamental modular region with 
the Gribov horizon). At least, some of these field configura- 
tions, as e.g. the center vortices, are topologically non-trivial. 
If the Gribov-Zwanziger scenario is correct they provide con- 
finement. But they may also provide a non-vanishing topolog- 
ical susceptibility and thus the C/yi(l) anomaly, especially the 
77' mass. When using, instead of a lattice Monte-Carlo calcu- 
lation, a functional approach (as Dyson-Schwinger or Renor- 
malization Group Equations) topological effects are encoded in 
the infrared behaviour of the Green functions. In Landau gauge 
it is the chain of infrared divergences of the ghost propaga- 
tor, the three-gluon and the quark-gluon vertex which seems to 
be responsible for both, confinement and the C/4(l) anomaly. 
This is, however, only a reflection of the fact that certain types 
of gluon field configurations cause both, confinement and the 
C/4(l) anomaly. 

If one is willing to accept this scenario a further conclusion 
is evident: the C/4(l)-relevant, topologically non-trivial gluon 
field configurations are intimately related to confinement, and 
thus are very likely not the instantons discovered already more 
than thirty years ago IfTSl . In addition, we know at least one 
type of gluon field configuration satisfying the requirements of 
this scenario, namely center vortices. However, at this point the 
relevance of other configurations to confinement and the Ua ( 1 ) 
anomaly cannot be excluded. 



7 Summary 

We investigated the Kogut-Susskind mechanism for the resolu- 
tion of the problem of QCD, employing known results 
for the ghost and gluon propagators from Landau gauge stud- 
ies of the Yang-Mills sector, together with recent results for 



the quark-gluon vertex. Taking the qualitative features of this 
study, we modelled the vertex-dressing such that an infrared 
soft-singularity was present, whilst at the same time match- 
ing meson observables for the vr and the p. The presence of 
a collinear singularity |8|, dominant in the gluon momentum, 
proves essential to the generation of a non-zero topological 
susceptibility in the chiral limit, as required by the Witten- 
Veneziano mechanism. Our study illustrates the effects of topo- 
logically non-trivial field configurations on the level of corre- 
lation functions and therefore provides a qualitatively feasible 
mechanism for the resolution of the mass of the 77' in a func- 
tional approach to QCD. 
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